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Abstract 


This is the second semi-annual report for NAS5-97039. Work in this project is 
related to NAS 1-20666, also funded by NASA ACMAP. The overall objective of this 
project is to improve the understanding of coupling processes between atmospheric 
chemistry and climate. Model predictions of the future distributions of trace gases in the 
atmosphere constitute an important component of the input necessary for quantitative 
assessments of global change. We will concentrate on the changes in ozone and 
stratospheric sulfate aerosol, with emphasis on how ozone in the lower stratosphere would 
respond to natural or anthropogenic changes. 

The key modeling tools for this work are the AER two-dimensional chemistry- 
transport model, the AER two-dimensional stratospheric sulfate model, and the AER three- 
wave interactive model with full chemistry. We will continue developing our three-wave 
model so that we can help NASA determine the strength and weakness of the next 
generation assessment models. 
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(I) Two-dimensional Chemistry-Transport Model 

(I I) Model Enhancements 

The following model enhancements were implemented in the past 6 months, 
primarily for the purpose of performing the WMO/UNEP ozone trends calculations. (1) 
We downloaded from the Goddard DAAC reanalyzed temperature data from NCEP for the 
1987-1995 period. This complements the 1979-1986 data we had obtained in 1996. The 
data was processed into monthly average zonal means and monthly temperature probability 
statistics on the AER 2-D CTM grid. We have used the data both to generate climatological 
temperature and temperature distributions for the entire 1979-1995 period, and for 
calculation of historical ozone trends employing the actual temperature for each year 

between 1979 and 1995. (2) The reaction CIO+OH— >HCl+0 2 was added to the 2-D CTM. 

Though it has been know for some time that use of this reaction would improve agreement 
between models and measurements near 40 km, a definitive measurement of the rate 
constant was not made until 1997 when Lipson et al. (1997) measured the branching ratio 
for the CIO+OD reaction and found 5-6% branching to DC1. (3) Five additional CFC and 
HCFC species were added to the model in order to account for future CFC substitutes in 
our ozone trends calculation. These were CFC-1 14, CFC-1 15, HCFC-141b HCFC- 
142b, and HCFC- 123. 


(1.2) Studies and Publications related to the 2-D CTM 

(1.2.1) ODPofCH TSr 

We submitted a paper to the Journal of Geophysical Research on the ozone 
depletion potential (ODP) of CH 3 Br. The ODP of CH 3 Br is determined by combining the 
model calculated ODP/BLP ratio (where BLP is the bromine loading potential) for CH 3 Br 
and its atmospheric lifetime. This paper examines how changes in several key kinetic data 
affect the ODP/BLP ratio. The key reactions highlighted in this study include the reaction 
of BrO+H0 2 , the absorption cross-section of HOBr, the absorption cross-section and the 
photolysis products of BrONO a , and the heterogeneous conversion of BrONO, to HOBr 
and HN0 3 on aerosol particles. By combining the calculated ODP/BLP ratio with the latest 
estimate of 0.7 year for the atmospheric lifetime of CH,Br, the likely value of ODP for 
CH,Br is 0.39. The model calculated concentration of HBr (~ 0.3 pptv) in the lower 
stiatosphere is substantially smaller than the reported measured value of about 1.5 pptv. 

The model can reproduce the measured value of 1.5 pptv if one assumes a yield for HBr of 
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1.3% from the reaction of BrO+OH, or a yield of 0.6% from the reaction of BrO + H0 2 . 
Our calculations show that the effect of these assumed rates on the model calculated 
ODP/BLP ratio is minimal: practically no impact for the assumed BrO+OH yield and 10% 
smaller for the BrO+H0 2 case. 

A copy of the manuscript is included as an appendix. 

(1.2.2). Adjusting Ky yin 2-D CTM according to observations 
a) Using exchanging rates to adjust in subtropics 

The subtropical barrier reduces the mixing between the tropics and the mid-latitudes 
in the lower stratosphere, and has significant effects on tracer distribution in the 
stratosphere. It is important for 2-D CTMs to have a realistic exchange rate between the 
tropics and mid-latitudes to simulate the effects of the subtropical barrier. The exchange - 
rate is determined by the eddy diffusion coefficient K yy in the model. In an earlier study, 

we introduced a tropical pipe by reducing the K yy in the tropics from 0.3 to 0.03 (10 6 

m 2 /sec), (Weisenstein, et al. 1996). Recently, in situ measurements of chemical species 
with a wide range of local lifetimes have been used to quantify the air exchange rate 
between the tropics and mid-latitudes in the lower stratosphere from the tropopause to 
around 21 km. It is found that the mid-latitude air is entrained into the tropical lower 
stratosphere with a replacement time scale of 10-18 months (Minschwaner et al., 1996; 

Volk et al., 1996). Meanwhile, Schoeberl et al. (1997) estimate the mixing time between 
the tropics and mid-latitudes in the 20-28 km region to be at least 18 months using the QBO 
signals in N 2 0/CH 4 ratio and tropical winds from UARS measurements. These results 
demonstrate that the tropical barrier allows a moderate penetration. 

The horizontal diffusion fluxes and the mixing ratio of CCI4, CF2CI2, CFCI3, 
CH3CCI3, CH4, HNO3, N2O, NOy, and O3 calculated in the AER 2 -D model are used to 
derive the exchange rates between the tropical and the mid-latitudes in the lower 
stratosphere. They are compared with the exchange rates estimated from observations. 

Using the original model (global diffuser model) with a K yy of 0.3 x 10 6 m 2 /sec in the 
lower tropical stratosphere the model calculated exchange time scales are around 5-6 
months. The tropical pipe model results show a much slower exchange rate of 38-60 

months. Our studies showed that values for K yy of 0. 13 x 10* m 2 /sec from the tropopause 
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to 21 km and 0.07 x 10 6 m 2 /sec from 21 km to 35 km in the model (leaky pipe model) will 
generate an exchange time scale of around 14 months below 20 km and 18 months in the 
20-30 km region, in good agreement with the estimates from observations. 

The model calculated lifetimes for species like N 2 0 and CFCs in the leaky pipe 

model are reduced by about 10% compared to the global diffuser model. Our model results 

show that the calculated ozone response to HSCT engine emissions is also smaller in the 
leaky pipe model. 

The paper that summarizes the results of this study was sent to Journal of 
Geophysical Research for publication. Two reviewers have made some suggestions and 
we are in the process of revision. The revised version was sent back to Journal of 
Geophysical Research in February 1998. 


b) Using latitudinal gradient of N 2 0 and NO/O3 to adjust near polar vortex 

It is well know that the isolation of the polar vortex is one of the preconditions for 
the Antarctic ozone hole to be developed (WMO 1988). The mixing ratio of long-lived 
tracers like CH4 and N20 usually are quite different inside and outside of the vortex due to 
the isolation of the vortex and the descending motion of the air inside the vortex during the 
polar winter night. The potential vorticity (PV), which is a good dynamic tracer, has been 
used to define the boundary of the vortex where the PV has large horizontal gradients 
(Manning et al., 1995). The gradient of the long-lived tracers across the vortex is usually 
underestimated in 2-D model simulations. There are two dynamical parameters in AER 2-D 
CTM which affect the horizontal gradient of long-lived species around the vortex, Ky y and 
the circulation. We found that the value of currently used in AER 2-D CTM during the 
winter high latitudes is probably too large to maintain sharp gradients for long-lived species 
near the polar vortex. We put a polar barrier in the model by reducing Kyy from 10 10 
cm 2 /sec to 10 8 cm 2 /sec between 50° to 70° during the winter. It is found that the horizontal 
gradient of N 2 0 and the ratio of NOy/0 3 near the vortex is quite sensitive to the change of 
Kyy. Fig. 1 shows the model calculated annually averaged ratio of NO y /0 3 at 64 mb with 
and without the polar barrier along with the observations. The polar barrier seems to have 
improved the model simulation. The generalized coordinate method, reported at the M&MII 
workshop (1997) by David Considine of GFSC, used to process UARS N 2 0 data 
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produces larger gradient of N 2 0 near the polar vortex than ordinary zonally average does. 
We will use this method to process UARS N 2 0, NO y and O, data and use their gradient 
near the polai vortex to fine-tune Ky y in the polar region. 


(II) 2-D Interactive Model 

(II. I) Improvements of the 3-wave Model 

The current version of the AER 3-wave interactive model has a vertical, resolution of 
3.5 km. In order to meet the future requirement for assessment studies, the resolution 
should be increased especially near the tropopause so that the exchange between the 
troposphere and the stratosphere can be better simulated. We decided to increases the 
vertical resolution of the model by three times to match the fine resolution of 1.2 km used 
in the AER 2-D CTM. The model with finer resolution is running on a SGI machine. The 
results show some differences in the temperature fields between two versions of the model 
with different vertical resolution. The causes of the differences could be the scale 
dependent parameters in the dynamics module such as the sub-grid fourth order diffusion 

coefficients. We are adjusting these parameters to improve the model with finer vertical 
resolution. 


The quasi-biennial oscillation (QBO) is the primary mode of variability in the 
tropical lower stratosphere and has been well documented (e. g. Reed et al., 1961). The 
QBO signals in the temperature have also been observed in the same area (e. g. Reed, 
1964). Associated with the tropical temperature QBO anomalies is a meridianal circulation 
m the tropics and subtropics, which should have a direct effect on the tracer distribution in 
the lower tropic stratosphere. QBO signals in total ozone have been observed for over 
three decades (Funk and Gamham, 1962). The distributions of N 2 0 and CH 4 observed by 
UARS/HALOE also show QBO signals (Schoeberl, et al., 1997). In order to improve our 
model’s ability to analyze the observed data or to use the model to assess atmospheric 
effects of pollutant, QBO should be included in the model. It is difficult for our 2-D 
interactive model to generate tropical waves to provide correct momentum sources for QBO 
due to the coarse vertical resolution and the simplified tropospheric dynamics. So instead, 
the zonal-mean QBO winds, averaged from NCEP data, are imposed on the model by 
adding a relaxation term in the zonal-mean zonal momentum equation. The time scale of 
lelaxation is about 100 days. The model generates the right QBO in the zonal winds but 
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only a weak QBO signal in the temperature. Accordingly, the QBO signals in ozone and 
other tracers are small. We are analyzing the model outputs to diagnose the causes of the 
weak QBO signals in temperature. 

(II. 2) Coupling of the 3-wave Model to the 2-D 19-Zone Climate Model 

In order to examine the total effects of solar variation and stratospheric aerosols, 
including their indirect effects from changes in ozone, we are upgrading the 2-D climate- 
chemistry model (CCM). One of the refinements is to replace the 2-D dynamic-chemistry- 
transport module of the 2-D CCM with the 3-wave model (Shia et al, 1997), which 
resolves the three longest planetary waves and calculates the eddy diffusion coefficient for 
chemical tracers from the resolved planetary waves. Thus, the model includes more 
interactions between chemistry, radiation and dynamics. In addition, the 3-wave model has 
updated heterogeneous chemistry. 


The new model is running well on a multi-processor SGI machine when the data 
exchange between the two modules, the 3-wave model and the climate model, is turned off. 
The next step is to establish full interactions between these two modules 

(III) Assessment Related Activities 
(III. 1) Models and Measurements Workshop (II) 

Our participation in the Models and Measurements Workshop has involved several 
of us at AER in different ways. Malcolm Ko contributed to the definition of one of the 
experiments. Courtney Scott is responsible for the analysis of the source gases 
experiment. The group, as a whole, has produced model results for the UADP archive. 
Malcolm Ko attended the first meeting in August and both Malcolm and Courtney attended 
the second meeting in November, 1997. At the meeting, Courtney made a presentation of 
the results from the source gases experiment. A write-up of the experiment was also 
submitted to UADP in January, 1998. 

Malcolm Ko will also help co-ordinate the preparation of the final report. 

(III.2) Ozone Assessment Report 
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Malcolm Ko is a co-author for Chapter 1 of the WMO Ozone Assessment Report 
with lead authors for the chapter, Ronald Prinn and Rudolf Zander. Chapter 1 deals with 
long-lived trace gases. Malcolm Ko worked with Michael Volk to prepare a section on the 
atmospheric lifetimes of the gases as can be determined from model calculations and 
observations. The model calculated lifetimes are taken from the Models and Measurements 
Workshop (II) calculations. The results from observations are taken from Volk et al. 

(JGR, 102, 25543, 1997) which uses historical records of tropospheric measurement to 
correct the observed correlation in the lower stratosphere to obtain steady state lifetimes. 
One result of the study is that there are indications that the lifetime of CFC-1 1 may be 
shorter than the currently accepted value of 50 years. A clearer picture should emerge 
when the uncertainty in the lifetime can be reduced. 

Jose Rodriguez is lead author for Chapter 2 along with Michael Kurylo. Jose’s 
effoit is funded by this project. Chapter 2 deals with source gases that are removed by OH 
in the troposphere and a large part of the chapter was devoted to methyl bromine (CH 3 Br), 
including its ozone depletion potential 

We also performed some chlorine loading calculations and communicated our 
results to Guus Velder who is one of the lead authors of Chapter 1 1 . Finally, results from 
a 1970-2050 trend calculation with the AER 2-D CTM were submitted to Chapter 12 by the 
January 1998 deadline. This calculation was performed with observed monthly mean 
temperatures and temperature probability distributions between 1979 and 1995. 


(IV) Chlorine loading and stratospheric cooling 

Recent winter/spring observations in the Arctic [Manney et al., 1996,1997; Muller 
et al., 1997; Rex et al., 1997] reveal further reductions of the ozone content despite nearly 
complete disappearance of Pinatubo aerosol from the stratosphere [Thomason et al., 1997] 
and essentially constant levels of chlorine and bromine loading [WMO, 1995]. Existing 
photochemical models of the stratosphere with heterogeneous chemistry predict little or no 
further stratospheric ozone depletion, unless an alternative mechanism is invoked. One 
such mechanism could be the cooling of the lower stratosphere [Oort and Liu, 1993; 

Randel and Cobb, 1994; Ramaswamy et al., 1996] leading to triggering of polar 
stratospheric cloud (PSC) formation. As a result, chlorine activation by enhanced PSCs 
provides further ozone decline during the winter/spring period. Here, we present sensitivity 
study results of the AER box model for typical conditions in the lower stratosphere at 70°N 
during the December 1 - April 1 period for an idealized parcel with different assumed 
stratospheric coolings and chlorine loadings. Our calculations show that the imposed 
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stratospheric cooling could further deplete ozone and retard its expected recovery even with 
the projected chlorine loading decrease. We show that, for the present conditions, a 1 K 
cooling could provide the same local ozone depletion as an increase of chlorine by 0.40- 
0.95 ppbv for the scenarios considered. Thus, sustained stratospheric cooling could further 
l educe Arctic ozone content and delay the anticipated ozone recovery in the Northern 
hemisphere even with the realization of the Montreal Protocol and its Amendments. 


(V) The POLARIS Campaign 

Malcolm Ko was in Fairbanks, Alaska for the deployment of the Polaris Campaign 
in July. Jose Rodriguez joined the September deployment, using the opportunity to meet 
with Michael Kuryio to work on Chapter 2 of the WMO Ozone Assessment Report. 

One focus of the analysis continues to be the NO s /NO y ratio. The observed 
NOj/NOy is larger than that calculated using a constrained photochemical model. R -S 
Gao (NOAA) and Ross Salawitch (JPL) have tried to use a rate constant for OH + NO, that 
is 50% slower than the JPL recommendation. We used the AER 2-D CTM to perform 
some calculations to see its effect on model calculated ozone trends. 

A copy of the manuscript is included as an appendix. 

(VI) Anticipated Activities in the Next Six Months 

Activities for the next six months will be dictated by the various assessment related 
activities of both the Models and Measurements Workshop, and of the WMO report. We 
hope to find time to revise some of our manuscripts for publications. 
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Where p is molecular weight in grams, and the mass of the atmosphere is 4.8x102' gm. The 
BLP IS the ratio of the steady state concentration of bromine released by compound x for a 

given emission rate to the steady-state concentration of chlorine released by CFC-1 1 for an 
equal (by mass) emission rate: 
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The ODP is defined as the ratio of the ozone depletion produced by a unit-mass 
emission of compound X to the ozone depletion produced by a unit-mass emission of CFC-1 1 
(the traditional reference gas). Numerical models that incorporate the pertinent chemistry and 
trace gas transport of the atmosphere are used to compute the ODP. For such modeling 
exercises, the steady state ODP for compound X can be given by (Fisher et al„ 1990); 
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cotTesponds to the spatial and seasonal average of the i i 

hrn . sonaj average of the local ozone removal efficiency of the 

7* re "*“ d * x — » - , *. , h „ 

, I' 0. if (wo species rele,,, bremi.e radio,,, « 

” - BEF valves will be drffeten, „ , wo 

d, ffc 7“ ” f " ““ — - ealculare *** 

affected 77“ “““ ^ '« “» »- specie, « 

So| ' “ r " 0, * , “ «“ “ 10 mode,. P„„oa, analyses 

f 3l '’ I992) h3Ve Sh ° Wn that th£ m0de ' ca,culated ODP/BLP ratio depends on the 

I “7 - — » lifetimes dne to tropclpbede « 

mix . S W ° Uld n0t affeCt thC ° DmhP ratio as <°"g « they are not so fast that the 

,nd - ° f * — ~ *— — - - - — ^ „,e o, ,;r 

propetty, one can define tbe ODP and BLP ,b, „co„mended ^ ^ __ 

Hows. The model is used to calculate atmospheric lifetime (T model ) and (ODP) modeI 

*"* «*“ *“ stmtospheric pr.oesses am 

accounted for. A value for fRT , 

r (t$LP) mode j IS then defined as 


(BLP) mo del = 


p? 

l 

[ T modeJ 

J. 

[S?]c 

FC-ll [ 

r model]c 
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and the model calculated BEF 


x ls given by 


BEF X = 

( BL n modcl 


The ODP value based on Ihe 


recommended lifetime is defined 


as 


(ODP)^rec - (BLP) Xrcc ( ^2^nxo6ei 

<< BLP )mo&el (BLP)TrecB ^ 


where (BLP) Trec j s defined 


as 


f-i'T] r -I 

(BLP) Trec = -AjL7l_._firecJ 

f~jr]cFC-M [tree] 


CFC-ll 


and 




Table I summarizes how the calrnR^ 1 

e calcu Dted values of ODP for OR p r u , 

1 990. These valine « , , C H 3 Br has changed since 

nese values are calculated by the AFR mnH.i • 

at the time. Between 1990 and 1997, the BEF for CH Br hi inc^ ^ ^ 

decreased. On balance, the ODP has decreased ’ mCreaSed Wh ‘' £ “ S 


from 0.99 to 0.39. 



Atmospheric Lifetime of Methyl Bromide 
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The net atmospheric lifetime (x lol ) for methyl bromide i 


is given by: 


'gas 




'land 


above expression, ^denotes the lifetime against atmospheric removal by gas-phase •' 

“ Ph “° ,ys “ ° rCH)B '' - “ — * - « «— — » « dea „ 

ntegrated loss rate (cnrV.) over the who)e atmosphere. The lifettme aga.nst ocean 

uptake, x ocean , is defined as the methyl bromide 


atmospheric burden divided by the flux i 


into the 


ocean. A similar definition holds for the lifetime 


against uptake by land surfaces, x !and . 


inth S P ° t0 ySIS and react i°ns of c H 3 B r with O('D), 0 ( 3 p )t NO 3 , and Cl take place 

c a mosphere, gas-phase removal of methyl bromide is likely to be dominated by its 
eactron with tropospheric OH. Based on the latest recommended rate constants of th ' 

ofCH 3 BrwithOH,andofCH 3 CCl 3 withOHto m • „ *e react, on 

ALE/GAP P j , D ' Wlth the “ ost recent analysis of 

CjAGE data fPrinn et pi iqq^n 

. " parttal lifetime due to reaction with OH in the 

troposphere for CH Br of ? i + n < - ne 

due ,0 stratospheric removal, the estimated x gu is 1.7 years. 


cT *“ - - - — « - ** 

C Hl B„ Tdls . . «* ,0 „ propofa , whicl duc 

““ y <W “' ly 



10 9 S ' ,) - haVC OCCan Up,ake lifeti " 1CS order of hundreds of years (Wine and Chameides. 
1990). Butler (1994) described a method to calculate the (partial) atmospheric lifetime of 

CH 3 Br due to removal by the ocean ( W in terms of the surface mixed layer depth of the 

ocean, sea-surface temperature, wmd speed, thermocline diffustvity, solubility of CH 3 Br in sea 

water, and degradation rate coefficient in sea water. The Butler (1994) estimate was carried out 

ustng a two-box model so that weighted averages of the parameters were used in the 

calculations. The calculated best estimate for t 77 

Wan Was 3.7 years, with a possible range from 

1.4 year to 14 years. These values were adopted in the WMO 1994 report. The same 

methodology was used in Won and Butler (1996) but with the parameters calculated 

dually for each 2 by 2 gnd over the oceans. Thus, they were able to take into account 

the co-variations of the parameters before summtng the effects to obtain the global partial 

hfetime. W.th this method, they obtained a best estimate for x ocean of 2.7 years, with a range 

from 2.4 years to 6.5 years. The latest update (Yvon-Lewis and Butler, 1997) includes 

biological removal and lowers the lifetime to 1.8 years, with a range from 1.1 years to 3.9 
years. 


Measurements of deposition velocity of CH 3 Br over land are not available. Based on 

the rat6S meaSUfed f ° r fiVe typeS of soi >. Shorter et al. (1995) computed the partial 
I.fetime of CH 3 B, The value is based on the sum of the annual uptake fluxes (g yH) by the 

five soil types, each of which is estimated as the product of the uptake rate for the soil type (g 

01 2 d ~‘ ) ’ the 3rea ° f the S0U type ° n the « lobe ( m2 >- ‘he number of active days in a year (d 
VH). The five soil types were assigned two uptake rates, corresponding to deposition 

ve.octt.es of 0.12 cm r « and 0.02 cm r«. The value obtained duoto soil removal is 3.4 years. 


Combining the three partial lifetimes 


provides a best estimate of 0.7 


years. 
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Partitioning of the Bromine Radicals 
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We used the AER 2-D model (We<sesntein et ah, 1996) with the standard JPL reaction 

r : COnStamS f ° r ,he gaS ' PhaSe r — parameterization for heterogeneous reactions on 
“ ate aeroso,. Foi, owing the nrethod used in WMO (.992), the OOP ls calcuiated based 
ackground chlonne and bromine concentration of 3.5 ppbv and 19 pptv, respectively. 


on a 


jpj Q7 6 COnCemra ' ,0nS ° f thg bromine radicals « 38»N for spring conditions ca.culated using 
recommended rate constants are shown in figure .. Note that at that latitude, BrONO, 

rO are the most abundant species. The concentrations of HOBr and HBr are less than 1 
pptv. Figures 2a and 2b show the computed [BrO ]/[Br 


by Avalonne et al. (1995). It is evident that 


larger than those calculated 


r y ] ratio along with observations reported 
the BrO concentrations calculated using JPL -97 are 


using JPL-94. This accounts for the larger BEF in Table 1. 


The reported measured concentration of 1.5 pptv for HBr between 25 and 35 km 

(Carlotti et al., 1996) is much larger than the model computed value of about 0.3 pptv. In 

section 6, we will discuss ways of increasing the calculated concentration of HBr in the model 
and how they affect BEF. 


5. Sensitivity of BEF to Rate Data 


The following reactions represent the primary bromine-driven 
cycles in the stratosphere: 


catalytic ozone removal 


CIO + BrO + hv Br + Cl + 0 2 

Br + 0 3 -> BrO + 0 2 
Cl+ 0 3 -> CIO + o 2 


(D 



(I) 
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and 


CIO+ BrO + hv -> Br + Cl + 0 2 
Bi + BrO + 0 2 

ci+ o 3 -> cio + o 2 

2 O 3 -> 30 2 . 


BrO + H0 2 — > HOBr + 0 2 
HOBr + hv — > OH + Br 
Br + O 3 — » BrO + 0 2 
OH + O 3 -» H0 2 + 0 2 , 


2 O 3 — > 30 2 . 

BrO + O > Br + 0 2 1 
Br + O 3 — > BrO + 0 2 J 

O +O 3 20 2 . 

2(Br + O 3 — > BrO + 0 2 ) 1 
BrO + BrO 2Br + 0 2 J 

20 3 -> 30 2 . 


BrN 03 + hv Br + NO 3 
NO 3 + hv NO + 0 2 
Br + O 3 — > BrO + 0 2 ^ 

NO + O 3 -> N0 2 + 0 2 
BrO + N0 2 + M -» BrN0 3 + M J 


20 3 30 2 . 


(II) 


(HI) 


(IV) 


(V) 


Current rate constants ind.cate that, in the present day atmosphere, the first three cycles 

contribute about an equal amount to ozone depletion in the lower stratosphere, with the last two 
about an order of magnitude smaller. 
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To ilJusuoie bow changes in * BEF, « ihe key 

““ JPL '’ 4 ” d " ” 7 T *“" 1 CI “"S“ » "»»"• <* ' > «™„gh <R S > (see Tahie 

> (he e„«,ey e,,h„ to „se (he, « (he „„ hnh.ing 

‘ ‘ Cy °"' “ " d '“' y ^ *» par,i(ioning he»« n BrO and ,he W, 

esplore ,h« d epe„d.„ee by pe„„™„g , ses« of c„c„„,io.s ,h, BEF in each case, 

ihe results are summarized in Table 3. 


Reaction R1 in Table 2 has received much 


attention. A key revision in JPL-92 that 


an estimated temperature dependence, the rate recommended in JPL-92 based on room 

temperature value by Poulet et al. (1992) ts about a factor of ,0 faster at stratospheric 

temperatures than that recommended by JPL (1990) i 

oy JrL, (1WU). Reaction R1 represents the rate limiting 

steps in cycle (II). With the rate constant from JPL-92 (a ten-fold 
90), the BEF increases from 37 to 49. 


increase compared to JPL- 


In the JPL-94 calculation, we included the absorption cross-section for HOBr as 
reported by Orlando and Burkholder (1995), which is the first measured absorption cross- 
section. P„or to that, the much larger cross-section of HOC1 was used as a proxy, resulting in 
computed HOBr concentration of less than 0.1 PP tv at 16 km and 0.8 pptv at 30 km. With the 
r ando and Burkholder cross-section, the concentration of HOBr is around 1 pptv at 16 km 
- 3 pptv around 30 km (see the curve labeled JPL-94 in Figure 3), comparable to those of 
and BrONO,. Inclus.on of this cross-section lowers the BEF from 49 to 44. 


Several revisions in JPL-97 identified in Tahlo 7 „ i ,. 

uentuied in lable 2 also have impacts on the model 

calculated HOBr. A smaller rate constant for R1 results in smaller HOBr (see Case B in Figure 
The „ oon> „ BfONO _ to BOBr and BiO. As a of 

(«,(». b«b „OB, -d „ l„g„ (« Case D in F„« 3). Tb, effect of R , 



"* “ “ 3 * C “ E F” « perfoTOd c „ „ w „„ we 

pdated all the JPL-97 rale data except the absorption cross-sections. The value obtained 

1,1 the l0Wer StratOS P here is Slm " ar to Case E. The difference between Case E and Case F 
around 30 km is the result of the slower rate for R8 (O + HOBr) in the JPL-97 

recommendation. Figure 4 shows the absorption cross-section of HOBr from Orlando and 

Burkholder (1995) and the JPL-97 recommendatton. As noted in JPL-97, the recommended 

cross-section is influenced by the measurements reported by Rattigan et al. (1996). With the 

larger cross-section between 300 nm and 400 nm and the extension of the cross-section beyond 

400 nm, the photolys.s rate of HOBr is three times larger in the lower stratosphere leading to a 
smaller concentration of HOBr (see cases F and JPL-97 in Figure 3). 


We will now discuss the effects of the rate revisions on the model calculated BEF for 
CH 3 B, The Orlando and Burkholder cross-sect, on decreases the sensitivity of the model 
ted BEF to changes m the reaction rate constant for reaction (1). With HOBr being a 
major reservoir, any increase (decrease) in k, is compensated by decrease (increase) in (BiO] 
so that the change in the rate limiting step k , [BrO][H0 2 ] wi„ be smaller than the comesponding 
change in k,. This is illustrated by the results for case B in Table 3. Despite a factor of two 
decrease in the rate constant for RJ, the OOP for case B is practically the same as that for case 
A. Revision of the BrONO, cross-section involves both the magnitude and the product yield 
The measurement reputed by Burkholder et al. (1995) extends the measured values to beyond 
390 nm. With the new cross-section, the photolysis rate is a factor of 1 .25 larger in the lower 

stratosphere. For BrONO, photolysis fR2l on Ivcb r 

y Br + N03 eac * s t0 ozone destruction as the 

°B^N02 branch is a null cycle. Previously, it has been assumed that 4> with * 

=0. The new values from N.ckolaisen and Sander (private communication) indicate that O 

Br 

= 0.29 with <E> —f>7i "n 

bk)+no 2 U/L Ihus the new branching should decrease BEF. However, 


BrO + N02 


+N03 
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diagnostics reveals that the cycle (V) contribution is small compared to the other cyc.es around 
20 km. Thus, the effect of the change in branching on BEF is small as illustrated by Case C. 
However, the expected effect from the extension of the absorption cross-section to longer 
wavelength is to suppress BrONO, and increase BrO and BEF (see Case F and Case H). 


Reaction (R3) is a major source of HO x in the lower 


stratosphere via production of 


HOBr and its subsequent photolysis. The reaction also leads to a higher concentration of BrO 

and smaller concentration of NO, The smaller NO.-provides the coupling to the C10 x and HO 

cycles. Cases D shows a 19% increase in the BEF. Our diagnostics reveal increased ' 

contributions from cycles (I) and (II) where a larger reseivoir of HOBr and BrO are produced 
via reaction (R3). 


The combined effects of Rl, R2b and R3 is given in case E. The resulted BEF can be 
compared with that for Case F, whrch consists of al, JPL - 97 updates except (he ___ 

e simrlarrity m BEF for the two cases shows that the effects of R5 through R8 on BEF are 
small. The final four cases in Table 3 illustrate the contributions of the revisions the cross- 
sections on BEF. The HOBr revision has the largest effects among the revisions. The larger 
cross-section leads to larger calculated concentrations of BrOand larger BEF. For BrONO the 

revision to the abso.pt, on cross-section again leads to a repartioning of the bromine radicall 

leading to higher concentrations of BrO. The effect from ih» „i, • , 

ne ettect from the change in branching is small. The 

cross-section revisions for CIONO and Rri-i 

° 2 and BrCl are each responsible for increasing the BEF by 

0.5 to bring the BEF to 58 in Case J. 
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The ODP and the BEF of CH 3 Br with HBr 


concentrations at 1.5 pptv 


The calculated concentration of HBr in Figure 1 is much 


smaller than values reported in 


"(e raw re. To „b,»„ dorUod »«„„„,<« of HBr similar ,o recenl w„ of 
HBr in the fowe, »,»„„« hy Carforri ,„ d Noll „ 

alternative mechanisms to produce HBr. 


Several studies (Larichev e, al„ 1995 ; Mellouk, et al.. !994) tndtcated that the product 
of R I ,s HOBr and 0 2 , although it is not possible to confinn or refute a small yield of HBr by 
Street measurements. Another posstbtlity is a small yield of HBr from the react.on of BrO with 
OH. Two recent publications presented model results to show how a small yield from BrO + 
OH (Chipperfield et a,„ 1997) or ytelds from either react, ons (Chartrand and McConnel 1998) 
could successfully stmulate the measured concentration, m two separate cases, we assume a 
1 -3% Yield of HBr from BrO + OH, or an HBr yield of 0.6% from BrO + HO, Annual 
averages of HBr between 25 and 35 Ion are calculated to be 1.3 pptv and 1 .5 pptv for the two 
pectively (see figure 5). In both cases, HBr is elevated by over 100% from Case A 
The concentrations of the remaining bromine species are decreased by up to 1 0%. The effect 
on BEF is small as evrdent from the small changes in BrO in figure 6. The calculated BEF 
values are 56 (for the BrO +OH case ) and 52 (for the BrO + Hn t 

„ , « a^ror the BrO + H0 2 ) case, compared to a value of 

calculated with no assumed HBr yield. 



Concluding Remarks 
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The calculated OOP for CH 3 Br depends not only on stratospheric bromine chemistry 
- also on the calculated atmospheric hfetune. Some of the processes which define the surface 
removal are not ye, wel, understood. In contrast to the vast amount of data on CO spectes 
o servations of BrO, species are few. Although observations of BrO in the lower stratosphere 
an he tentat.ve upper limits on HBr above 30 km place some constramts on bromine 
chemistry, ihey do „„ „ le model, which lissome , , m ,|| (5 . 10%) HBr ta|> , 

Laboratory amd.es o, , M (2) ^ ^ 

(Ota, Kilnde legion, wl„ help accelerate on, understanding ot t , ic 

chemistry.' 


bromine 
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3rN n, ■^zz'z'^r r pr °' ,i “ r ° r •** —«■» « 

were performed using JPL -97 r ,,„ 

fa 97 rate recommendations. 

F *“ 2 ‘ : C ° mp '“" •*»>« value, o, fc Br0m 

values derived from the AASE n me ^ C ° mpared with the 

— — - - » ta - *■ «. 

calculated value, of B,Y „ d The ra “ C ™ "»g model 

results are sorted by latitudes. 

F ‘gure 2b: Same as figure 2a 


except that the results 


are sorted by solarzenith angles. 


Fl gure 3: Model calculated dav f 

TaWe 3 description of ^ COnCentra ‘ ,0n * ^ ** -es. See 


Figure 4: Absorption cross-sections of HOR r 

JPL-97 recommendation. ^ ° ** Burkho,de r (1995) and the 


fyfode, calculated mixing ratio profiles for HBr 1 , 

"" "“”"»««■ ™e c„„o labeled Brf, „ OH " W4 “* ' PMl 

“ ‘ “* '«“«■ ™e curve Bl0 . m “ “>>'** 

the react, on. assumes an HBr yield of 0.6% f rom 


Figure 6: Same as figure 5 except for BrO. 
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I able 1: Values of ODP for CH3Br 


circa 

1992 

BEF 

36 

49 

key revision 

Tree 

(years) 

1992 

BrO+H0 2 10 times 
faster 

2.1 

2.1 

1994 

48 

7PLT94 ~~ 

1.3 

1994 

44 

HOBr cross-section 
from Orlando and 
Burkholder (1995) 

1.3 

1997 

58 

see table 2 

0.7 


comment 


removal by OH only lfO.73 
removal by OH only ~0 

include removal by "O' 
ocean 

include removal by "(X 
ocean 



include removal by 0. 

soil, and biological 
removal in ocean 
water 


Table 2: Summary of rate revisions for key bromine chemistry 


reactions in JPL-97 


RirBrO +HO, 


R2a: absorption 


section of BrONQ 2 , 


R2b: quantum yield 


R3: heterogeneous 
conversion of BrON0 2 
on sulfate 

R4: photolysis of 
HOBr, 


R5: BrO + O 


R6: BrO + BrO 


R7: BrO + OH 


R8: O + HOBr 


JPL-97 j jpl-94 


3.4x10-1 2 exp(540/T) 6. 2x 1 O' 1 2 exp(500/T) 


temperature dependent 
cross-section extended 


to longer wavelengths 


<£> =0 71 

BrO + N02 

<J> =0 29 

- BftNQl y 


temperature dependent no recommendation 


remark 


new rate is 0.65 times 


old rate @ 230K 


new J rate is 1.25 


times larger than the 
old one 


cyle V is 0.3 times 


smaller 


larger cross-section at 


long wavelengths. 


1.9x10-1 1 exp(230/T) 


1.5x10-1 2 exp(230/T) 


7.4xlO-U 


1.7x10"! lexp(260/T) 


4. Ox 1 0" 1 2 exp(- 1 90/T) 
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Table 3: Model calculated BEF using different input 
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Abstract 

Recent wmter/sprmg observations in the Arctic [Manney et a]., 1996,1997; Muller et a]., 

1997; Rex et al., 1997] reveal further reductions of the ozone content despite near complete dis- 
appearance of Pinatubo aerosol from the stratosphere [Thomason et al., 1997] and essentially 
constant levels of chlorine and bromine loading [WMO, 1995], Existing photochemical mod- 
els of the stratosphere with heterogeneous chemistry predict little or no further stratospheric 
ozone depletion, unless an alternative mechanism is invoked. One such mechanism could be the 
cooling of the lower stratosphere [Oort and Liu, 1993; Ilandel and Cobb, 1994; Ramaswamy 
et al., 1996] leading to triggering of polar stratospheric cloud (PSC) formation. As a result, 
chlorine activation by enhanced PSCs provides further ozone decline during the winter/spring 
period. Here, we present sensitivity study results of the AER box model for typical conditions 
in the lower stratosphere at 70°N during the December 1 - April 1 period for an idealized parcel 
with different assumed stratospheric coolings and chlorine loadings. Our calculations show that 
the imposed stratospheric cooling could further deplete ozone and retard its expected recovery 
even with the projected chlorine loading decrease. We show that, for the present conditions, a 
1 K cooling could provide the same local ozone depletion as an increase of chlorine by 0.40-0.95 
ppbv for the scenarios considered. Thus, sustained stratospheric cooling could further reduce 
Arctic ozone content and delay the anticipated ozone recovery in the Northern hemisphere even 
with the realization of the Montreal Protocol and its Amendments. 

Introduction 

It is now widely accepted that the spring ozone hole, first detected in Antarctica [Farman 
et al., 1985; Chubachi et al., 1985], can be explained primarily due to activation of chlorine- 
containing molecules of anthropogenic origin [Molina and Rowland, 1974] via heterogeneous 
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Initialization of the AER Box Model 

Photochemical box models have been applied successfully to analyze different aspects of the 
Antarctic ozone hole [e.g. Schoeberl et al., 1996; Danilin et ah, 1996; Kawa et al., 1997]. For 
this study, we use the AER photochemical box model with the heterogeneous chemistry package 
described in Danilin et al. [1996]. We investigate chemical effects in the lower stratosphere, 
since the ozone profile has its maximum there and its perturbation due to PSC processing is 
also strongest there, thus determining to a large extent the magnitude of total ozone column 
depletion. So our box model calculations in the lower stratosphere could serve as a proxy of 
the ozone column change. We simulate the behavior of an isolated air parcel at a latitude 
of 70°N subsiding with a rate of 1 km/month from. —24 km on December 1 to —20 km on 
April 1 [Rosenfield et ah, 1994]. The latitude of 70°N for our calculations is chosen because 
this latitude is typically located within the polar vortex according to analysis of zonal fields of 
potential vorticity [Manney et ah, 1994]. We constrain initial concentrations in the air parcel 
using the UARS data in a fashion similar to Carslaw et al [1997a]. Initial mixing ratios of 
Os ( 3 - 6 ppmv) , HC1 (1.7 ppmv), H 2 0 (5 ppmv), and CH 4 (1 ppmv) are chosen based on the 
HALOE measurements on December 1 at our starting point (70°N, 32 hPa). Initial values of 
C10N0 2 (1.2 ppbv) and HNO3 (11-5 ppbv) are taken from the CLAES measurements [Roche 
et al, 1994]. We assume 3 ppbv of Cl y (=C1 + CIO + HC1 + C10N0 2 + HOC1 + 2(C1 2 0 2 + 
Cl 2 )), 14 ppbv of NO y (= HN0 3 + NO + N0 2 + N0 3 + HN0 4 + HONO + 2N 2 0 5 + C10N0 2 ) 
and 20 pptv Br y (= Br + BrO + HOBr + BrON0 2 + BrCl + HBr) consistent with the CLAES 
measurements N 2 0 using the N 2 0-NO y [Loewenstein et al, 1993], N 2 0-Cl y [Woodbridge et al, 
1995], and N 2 0-Br y [R.Salawitch, personal communication] correlations. Initial partitioning 
among the rest of the chlorine, bromine, and nitrogen species is taken from the AER 2-D model 
[Weisenstein et al, 1996]. We also assume decrease of NO y in the air parcel from 14 ppbv on 
December 1 to 10.2 ppbv in late March to be consistent with MLS HN0 3 measurements [Santee 
et al, 1997]. Initial aerosol surface area density (SAD) is equal to 1.8 /rm 2 /cm 3 with assumed 
log-normal distribution with mean diameter of 0.155 ^m and width <7=1.49, consistent with the 
aircraft measurements for the background conditions used in Carslaw et al [1994]. The aerosol 
SAD is kept constant during the NAT runs and is variable during the STS runs due to uptake 
of HN0 3 and H 2 0. The shape of the temperature history chosen in this study (see Figure 1) 
mimics the evolution of the vortex-averaged minimum temperature in the lower stratosphere 
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[Zurek et ah, 1996] but is mcreased uniformly by 5 K to provide -30% ozone depletion at the 
pomt of interest during a model run from December 1 to April 1. This temperature shift serves 
as a ’tuning’ parameter. We also do not superimpose 6-10 day waves on our temperature history 
(like in [Carslaw et al., 1997a]) to avoid additional complications of the day-to-day variability 
of PSCs, which could mask the main point of our study. 


Type I PSC schemes 

One of the key uncertainties for our analysis is what PSC scheme should be adopted for the 
model runs. The issue of how type I PSCs are formed is still a subject for discussion [Carslaw 
et ah, 1994, 1997b; Tabazadeh et ah, 1994; Tolbert, 1994], However, the model analysis by 
Sessler et ah [1996] shows that, evaluating the chlorine activation and consequent ozone de- 
pletion, a whole spectrum of type I PSC formation possibilities is embraced by the nitric acid 
tnhydrade (HN0 3 »3H 2 0, NAT) and supercooled ternary solution (H 2 S0 4 /HN0 3 /H 2 0, STS) 
schemes. Our box model allows type I PSCs formation via NAT without supersaturation (with 
fixed drameter of 1 m) or via STS [Tabazadeh et ah, 1994], The comparative efficacy of chlo- 
nne activation in these schemes also depends on the assumed level of background sulfate aerosol 
and size of NAT particles. Figure 2 addresses differences between these two schemes from the 
pomt of view of SAD available for heterogeneous reactions in the NAT and STS schemes as 
well as the reaction probabilities and e-folding times of the principal chlorine activation reaction 
C10N0 2 + HC1 -4 Cl 2 + HN0 3 . A very steep increase of the reaction probability and surface 
area with decreasing temperature below 196 K drastically reduces the e-folding time of chlorine 
activation from -200 hours at 196 K to 1.5-3 hours at 191 K. For the particular parameters of 
this study, the STS scheme is more efficient below 191 K, while a steep increase of the the NAT 
SAD between 195 and 192 K makes the NAT scheme slightly more efficient than STS scheme in 
this temperature range. Figure 2 also shows that it is barely possible to use CIO measurements 
to distinguish the NAT and STS schemes for the air parcels exposed to low temperatures for 
a long time (days or so), since all chlorine is expected to be activated in both schemes. Addi- 
tionally, Figure 2 can demonstrate why the ozone hole has been first observed in the Southern 
(but not in the Northern hemisphere) in the 1980s. Colder temperature in the Antarctic lower 
stratosphere (by 7-8 K comparing with the Arctic stratosphere) and more stable polar vortex 
greatly facilitates chlorine activation, thus producing stronger ozone depletion for the same or 
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slightly lower level of Cl y there compared with that in the Arctic. 


Equivalence of stratospheric cooling and increase of chlorine loading 

To investigate the sensitivity of the ozone depletion to a wide range of chlorine loading and 
stratospheric cooling, we performed model runs for 2 to 3.5 ppbv Cl y and 0 to 5 K stratospheric 
cooling. Effects of stratospheric cooling are investigated by lowering the whole temperature 
pattern by 1 to 5 K while keeping its shape unchanged as shown in Figure 1. The change of 
ozone mixing ratio in the air parcel on April 1 as a percentage of its initial value of 3.6 ppmv 
is plotted as a function of chlorine loading and further stratospheric cooling in Figure 3. This 
figure summarizes our findings and underlines the widely accepted point that ozone content 
decreases with increasing level of Cl y (draw any horizontal line in Figure 3). However, this 
figure also clearly shows that possible further stratospheric cooling (due to C0 2 or H 2 0 [IPCC, 

1995; WMO, 1995]) has the effect equivalent to an increase of chlorine loading from the ozone 
balance point of view. 

In order to further elucidate results in Figure 3, we introduce “chlorine-cooling equivalent” 
(CCE) as a measure of the change in chlorine loading equivalent to a 1 K cooling for a given level of 
the ozone depletion in the air parcel under specified ambient conditions. Our results for the NAT 
scenario show, for example, that 30% decrease of the ozone mixing ratio in the air parcel could 
result from either 3.1 ppbv Cl y and T 6ase , or 2.15 ppbv Cl y and T 6ase -1K. In other words, a 1 K 
cooling produces the same effect in the local ozone content as an increase of Cl y by 0.95 ppbv 
(i.e. the cooling-chlorine equivalent is equal to - 0.95 ppbv(CI y )/K). For the STS scenario, the 
ozone sensitivity to cooling is smaller (CCE^0.4 ppbv(Cl y )/K for the 30% depletion isoline). 
The higher values of the CCE for the NAT runs can be explained by the very steep increase of the 
NAT SAD for the chosen present day case (see Figures 1 and 2). However, further cooling leads 
to the saturation of the NAT SAD due to complete conversion of HN0 3 from gas to solid phase, 
thus lowering the CCE values to several tenth of ppbv(Cl y )/K. The similar pattern is observed 
for the STS scenario, thus showing that the CCE concept developed here remains valid for the 
different PSC schemes. The absolute values of the CCE is a non-trivial function of the chosen 
ozone depletion, ambient atmospheric conditions, and chosen PSC scheme. The CCE is a local 
parameter and should not be extrapolated to the global scale. However, the concept of CCE, 
applied in more sophisticated global models, could be useful in regulatory discussions. One can 



also introduce a carbon dioxide - chlorine loading equivalent or ozone depletion potential (OOP) of 
C0 2 which could help to formulate a future regulatory ozone protection policy translating C0 2 
emissions to CFCs emissions after proper scaling using global model calculations. The WMO 
[1995] Report does not consider C0 2 as an ozone depleting species. 

Our calculations for Antarctic conditions (similar to [Danilin et al., 1996] and not shown 
here) show CCE values *0.15 ppbv(Cl,)/K there. This implies that possible future cooling 
of the Antarctic stratosphere would have a smaller effect on ozone depletion (which is close 
already to its saturation level [Hofmann et al, 1997]) than similar cooling in the Arctic. 


Discussion and outlook 

As a continuation of our sensitivity analysis, we present calculations of ozone behavior 
along the idealized air parcel considered in the previous sections using anticipated future level 
of chlorine loading and assumed stratospheric cooling. We follow the WMO scenarios for the 
Cl, (shown by solid black line in Figure 4) and CH 4 concentrations assuming a fixed level 
Of Br„ (-20 pptv). Figure 4 shows that, ignoring stratospheric cooling, one indeed should 
anticipate a steady recovery of ozone (green lines) in the lower stratosphere immediately after 
reduction of the Cl, loading peaking at the level of 3.08 ppbv at 24 km (or ~3.7 ppbv at the 
top of the stratosphere) in the year 2000. The red dashed line shows “the road not taken” 
[Prather et al., 1996], assuming the level of chlorine is fixed at 3 ppbv and stratospheric cooling 
with a rate of -0.15 K/year. However, even for the relatively modest constant temperature 
trend of -0.05 K/year [Randel et al., 1994; Ramaswamy et al., 1995], ozone remains below its 
1995-level for the next fifty years in the NAT scenario. For the STS scenario, ozone returns to 
the 1995-content near the year 2030. For the larger temperature trend of -0.1 K/year, further 
delay of the ozone recovery is expected with more than 40% ozone depletion, keeping the ozone 
content in the parcel below its 1995 level during next 50 years in the both scenarios. For the 
more severe cooling rate of -0.15 K/year, a larger ozone depletion is calculated (more than - 
40%) with a longer-lasting impact. For the considered scenarios, the maximum ozone depletion 
could be anticipated during the 2010-2020 period for all temperature trends and both PSC 
scenarios. Thus, our calculations show that the ozone depletion in the Arctic stratosphere may 
be deepened and its recovery delayed by several decades when stratospheric cooling is taken into 
account. Both our results and the 3-D climate model calculations of Shindell et al [1998] show 
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that the currently accepted anticipation of ozone healing in the Arctic following the elimination 
of CFGs emissions may be too optimistic. 

Our results show that the Arctic ozone is very sensitive to the ambient temperature in the 
lower stratosphere. We present estimates of only the chemical effects of stratospheric cooling. 
Stolarski [1997] discussed three possible scenarios of future development of the Arctic ozone, 
depending on what the main reason of the observed stratospheric cooling is: a) greenhouse gases, 
b) ozone, or c) just natural variability. Since the natural variability of the Arctic winter strato- 
sphere is large and relevant data sets are short, it is difficult to determine the cause of the 
observed cooling. Each cause could have important implications, a) Carbon dioxide provides 
the strongest cooling in the lower stratosphere among greenhouse gases (e.g. IPCC [1995], 
Clough and Iacono [1995]). If the cooling and the polar ozone hole propagate upwards in the 
near future, it could be considered a manifestation of greenhouse gas effects due to the vertical 
profiles of the C0 2 cooling rate [e.g. Clough and Iacono, 1995], In this case, the extent of 
polar stratospheric ozone depletion could become considerably larger and its recovery delayed 
by several decades even with the implementation of the Montreal Protocol. If future polar 
measurements and more sophisticated 3-D calculations support our findings, C0 2 emissions 
should be contained for better ozone protection. 

b) Ozone itself is a very important radiative species and there is a strong coupling between 
radiative, chemical, and dynamical processes in the Arctic stratosphere. In this case, more so- 
phisticated 3-D model calculations with all these feedbacks and a state-of-the-art stratospheric 
heterogeneous chemistry module [e.g. Austin et ah, 1992] are required to predict future cooling 
and evolution of the possible Arctic ozone hole. 

c) It is also possible that current stratospheric conditions represent a complicated mixture of 
natural variability and human impact. If this is true, further model prediction of the Arctic 
stiatospheric ozone recovery is an extremely complex issue. 

Since the current atmospheric conditions in the Arctic are very close to conditions allowing 
further steep ozone loss via an additional PSC formation, a very careful analysis of the con- 
sequences of current or future levels of emissions, which directly (e.g. aircraft H 2 0 and NO*) 
or indirectly (C0 2 ) facilitate the formation of PSCs or provide additional sites (aircraft sulfate 
particles) for heterogeneous chlorine activation, is required [Peter et ah, 1991; Weisenstein et 
ah, 1996; DelNegro et ah, 1997], Our previous study shows a possibility that severely depleted 
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ozone in the polar region could be diluted to the mid-latitudes after the final stratospheric 

warming [Sze et al., 1989], A more detailed study of the reasons for the lower stratospheric 

cooling with particular emphasis on anthropogenic contributions to this effect and consequent 

effects on ozone at polar and middle latitudes will require a use of a 3-D model with all neces- 
sary feedbacks. 
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Figure captions 

Figure 1 . The temperature history adopted for the box model runs at 70°N: solid black 
line show the present day case; black dashed lines depict the present day case lowered by 1, 2, 3, 
4, and 5 K; the almost horizontal dashed lines correspond to the NAT formation threshold and 

the ice frost point, respectively (their slope is due to a subsidence with the fixed H 2 0 mixing 
ratio of 5 ppmv). 

Figure 2. Top panel: Calculated surface area density of STS (solid red line) and NAT 
(solid blue line) particles in ^m 2 /cm 3 (left vertical axis) and reaction probability (right vertical 
axis) of the C10N0 2 + HC1 > Cl 2 + HN0 3 reaction on STS [Hanson and Ravishankara, 1994] 
(dashed red line) and NAT [Hanson and Ravishankara, 1993] (dashed blue line) as a function 
of temperature. Surface area of NAT was calculated assuming its unimodal distribution with 
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diameter of 1 and density 1.62 g/cm 3 . Below 205 K, the model of Tabazadeh et al [1994] 
is used to calculate swelling of aerosol particles due to additional uptake of H 2 0 and HN0 3 . 
Bottom panel: Timescale of the C10N0 2 + HC1 — > Cl 2 + HN0 3 reaction (defined as r 
= 4/(7xSurfaceAreaxSpeed CTCW 02, in hours) on the STS (red) and NAT (blue) shows an 
effectiveness of chlorine activation as a function of temperature. Black symbols show minimum 

temperature during the model runs from 195.25 K (present day case) to 190.25 K (present day 
case-5K). 

Figure 3. Isolines of ozone change in the air parcel on April 1 (in % from its initial value 
of 3.6 ppmv) as a function of chlorine loading and stratospheric cooling at 20 km and 70°N for 
the NAd (top panel) and STS (bottom panel) scenarios. These values correspond to constant 
Br y =20 pptv and the temperature evolution shown in Figure 1 (i.e. zero cooling correspond to 

T 6ase line and 5K cooling - to the T 6ase -5K line shown in Fig. I). 

Figure 4. Ozone depletion as function of year for the projected chlorine and methane 
loading and different rate of stratospheric cooling for the NAT (top panel) and STS (bottom 
panel) scenarios. Bromine loading is assumed to be constant at 20 pptv. The black line in 
the bottom panel shows the adopted decrease of chlorine loading at the starting point of the 
air parcel at 32 hPa. We considered scenarios when the present temperature pattern (Fig.l) 
remains. constant (green curves) and uniformly lowered at a rate of 0.05 K/year (blue curves), 
0.1 K/year (cyan curves), and 0.15 K/year (red curves). Ozone change in the parcel for the 
constant level of chlorine loading at 3 ppbv for the adopted cooling rate of 0.15 K/year is shown 
by the pink dashed lines. 
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